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The interaction problems of wave fields with resilient bounded bodies arise 
in the different fields of mechanics, hydro acoustics, geophysics and seismology. 

The analytical solutions of the problems of hydro acoustic waves scattering 
by resilient bodies are directly connected with the methods development of mod- 
eling wave fields interaction processes with underwater objects. Increasing the 
complexity of underwater objects geometry we get the increase of mathematical 
complexity of the pressure field determination problem solution. In this con- 
nection the problems are solved with simplified mathematical statement, which 
analysis helps to indicate the significant features of geometry and structure of 
modeling object. Up till now number of analytical and experimental researches 
were implemented according to the statement of the problem. These researches 
were directed to the conformities determination of influence of any construction 
or technologic property on structure parameters of reflected signals. 

The problems of scattering of acoustic pressure waves by resilient bodies 
with complex geometry have no analytical solution. 

Applying the direct numerical methods (e.g. finite difference method) to 
solve partial differential equations of these problems is not expedient even using 
powerful computers due to multi-dimensional type of the problem. For this 
reason high-frequency and low-frequency methods are often used to solve such 
problems. But when wavelength is comparable by order with the characteristic 
size of disperse body the frequency lies in moderate range. Exactly this range 
as shown in jl6| is highly informational. 

To solve the scattering problems in the moderate frequency range the limited 
integral equations method is applied. 

Applying Grin theorem let us reduce the solution of stationary problem of 
resilient body dispersion to the solution of combined singular equations for the 
system of potential displacements. The solution of the above determines the 
surface potentials, which are used to define the outside field. Thus, the num- 
ber of dimension variables is reduced and it becomes possible to use numerical 
methods. 
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In works ^| ^ the problem of sound diffraction on the finite cyhnder with 
mixed boundary conditions is solved by Grin's function method. In ^Sj the 
axis-symmetrical dispersed pressure of resilient disk is computed using finite 
element method and HelmHoltz integral equation. 

As is well known, there is a tight relation between the solution of scattering 
and emission problems. It could be defined by reciprocity theorem |19j . This re- 
lation could be used to transfer emission problem results to scattering problems 
and vise versa |2U| . 

While solving complex problems, for example three-dimensional problems of 
acoustic waves scattering on the resilient bodies or two-dimensional problems of 
scattering on the piece-wise-smooth resilient bodies and problems of scattering 
of acoustic waves on the thin-wall covers, troubles with definite satisfaction of 
the appropriate boundary conditions arise. 

The considered work is devoted to the research of scattered and emissed 
acoustic waves by finite cylindrical covers, resilient cylinders and round plates. 

The presence of ribs fracture line in bottoms interface and cylindrical covers 
on the face of resilient cylinder and on edge of round plates inputs the additional 
troubles of evaluating of such problems. 

During last years for the solution of space problems of scattering of acoustic 
waves numerical and integral equations methods are effectively applied. De- 
spite the achieved researches of applying these methods to three-dimensional 
problems, questions of efficiency increase of existing methods and their PC- 
realization are still actual. In particular, one of the hydro-acoustic problems is 
the evaluation of scattered and emitted pressure field by three-dimensional bod- 
ies of cylindrical form. But during the scattered pressure field exploration near 
bounded disperser with edges and also internal features of dissipater troubles 
arise, which are hardly solved without analytical research. That is why the im- 
portant problem of analytic-numerical methods development and development 
of software for acoustic space interaction with finite bodies of cylindrical form 
are actual. 

During the solution of problems of emission and scattering of acoustic waves 
by bounded bodies of cylindrical form we use the strategy based on applying 
of bounded integral equations method combined with series method and later 
using the improvement of convergence method, which consider particularities 
of desired functions. The proposed strategy is authorized by the fact that the 
series method is effective to solution of equations of resiliency theory, thin covers 
and plates, and the bounded integral equations method is effective to research of 
infinite acoustic environment vibrations. Using this methodology let us describe 
the behavior of bounded resilient bodies near salient line, explore the pressure 
field in near and far zone, obtain the solution of the problem with defined 
precision in wide frequency range. 

In the first chapter basic linear equations of motion and initial relations are 
considered, statement and solution of problem of scattering of fiat acoustic wave 
by finite cylindrical cover bounded at the ends by resilient bottoms is done. 

The improvement of convergence of obtained infinite systems of linear alge- 
braic equations was held, using asymptotic properties of Fourier decomposition 
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coefficients of tlie desired functions while solving the equations of cover theory. 

The solutions of the problem of scattering of flat pressure wave by open-end 
free cylindrical cover, and also by tension cylindrical cover with defined inner 
tensions. In case of thin resilient cylindrical cover and round plates its motion 
is modeled by the linear Thimoshenko model. 

Today the problem of learning of dynamical processes of acoustic fields inter- 
action with thin-walled structures in liquids is of great interest. The theoretical 
exploration of these processes is strongly associated with building of effective 
analytic- numerical methods of mathematical physics equation solving, which 
are realized on today's PC. For example, while analyzing wave diffraction on 
solid bodies and covers of revolution with fixed ends in j^J finite-element and 
spline- functions methods based on function approximation theory were applied. 

We assume that scattering objects could easily move in boundless liquid 
environment or to be deformless. 

In this work we propose the analytical method of dynamical problems solu- 
tion of fixed oscillation of infinite acoustic environment interacted with resilient 
finite cylindrical bodies. 

In the nowadays technics problems of oscillation of finite cylindrical covers, 
filled with liquid, their sound waves emission at vibration in infinite acoustic 
environment are actual. 

While interaction of wave fields with resilient bounded bodies problems arise 
in different fields of technics, medicine, hydro-acoustics and geo-physics. 

The analytical solutions of problems of scattering of acoustic waves by 
resilient bodies are directly related with developing methods of modelling of 
interaction processes of wave fields with underwater objects. 

Let in the infinite ideal compressible rigid media the cylindrical cover of 
finite length bounded at the ends by bottoms is placed. Cover is filled with 
ideal compressible fluid inside with density po and sound velocity cq. 

The considered cover is referred to cylindrical system of coordinates r' , z' ,Q 
where axis oz' is aligned with axis of cylindrical cover and plane z' = is 
equidistant to ends (see Fig. 1). Flat harmonic in time pressure wave falls to 
the cover (see Fig. 1): 

(r, z, O, t) = Pq exp {—iuj (r cos 8 sin ip* + z cos (p* — t)) (1) 

In formula and below following symbolism is used: = "^5^ = 
nondimentional coordinates, t — -^:uj — are nondinientional time and 
angular frequency, c is sound velocity in the dissipating body environment, 
Rq is radius of middle surface of cylindrical cover, t, are time and angular 
frequency, if* is angle between wave direction and axis of the cylindrical cover, 
Po is constant with the dimension of pressure. Time factor exp (itor) is omitted 
below. 

The problem of definition of dissipated pressure field outside the cover P*^ (r, z, 9) 
and in filler P° (r, z, Q) and components of cover movement vector U and plates 
W is reduced to solution of system of differential equations. 
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Fig. 1 

Dissipated acoustic pressure (r, z,9) in unbounded acoustic medium is 
described by wave equation pP 

(A + (r, z, 0) = (2) 

where A is three-dimensional Laplace operator. 

Dissipated pressure in fluid which fills the cover P° (r, z, 6) satisfies the wave 
equation [2] 

(A + ^)p°(r,z,0)=O (3) 

where (31 = ^. 

The dynamics of thin resilient cylindrical cover is described by linear equa- 
tion of Timoshenko cover theory which takes into account rotary inertia and 
deformation of transversed shift 2 

LijUj = giSa; (Ly = = 1,2,3,4,5) (4) 

where 5ij is the Cronicler's symbol. 
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El, 1^01, Pi, hi is Young's modulus, Poisson's ratio, density and thickness of 
cover material, p is density of outside liquid medium, U2, U4 are turning angle of 
normal of cover's middle surface in planes rQ and rz appropriately; gr is radial 
component of outside force falls on middle surface area unit. 

Differential equations of Timoshenko-Mindlin which describe the transverse 
motion of round resilient plates with shifted and rotary inertia taken into ac- 
count and characterizing few first asymmetrical thickness mod oscillations have 
the following form |3j 

(A2 + ai Ao - 02) W, (r, 9) ^ (a* - Aq) (5) 



where 



{Ao + kl)^ir,e) = (6) 
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gz - normal component of outside force, which falls falls on middle surface 
area unit; -£2 , 1^02 , G2 , P2 , ^2 is Young's modulus, Poisson's ratio, shift ratio, 
density and plate thickness appropriately, Kt2 is numerical shift coefficient; $ 
is auxiliary function; C22 is velocity of transversed waves in plate. 

Let us define the radial and tangent components of motion vector of plates 
Wr, We which characterize the flat stress state as follows 



dip 19* 
dr r dr 



We 



Idif 9* _ 
r 86 9r ' 



(7) 



Here the scalar potential if and non-zero component of vector potential satisfy 
the wave equations 



Ao + ^j (^(r,0) = 0; (^Aq + 



(8) 



where 
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The solution of differential equations (0)-® must satisfy the following boundary 
conditions OIH] 

a) continuity translation on middle cover surface and plates 



d_ 
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b) continuity translation in junction place of cover and plates 
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c) continuity of turning angle (r = I, z = ±Z) 
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_Ro dr Ro 
d) continuity of forces (r = 1, z = /) 
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1 ^ i^oi / 1 dur, dui 
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dr r 



e) continuity of moments (r — 1, z — I) 
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Taking into account the symmetry relative to plane = we find the 
solution of equations ©-© in the Fourier series form: 

oc 00 

Ml (z, ^) = ^ ^ iyXnv COS f3i,z + u^^^ sii\(3yz) sinn^; 



oc 00 



U2 {z, ^) = X/ X/ ('"■2rw COS f3uZ + Mjny sin(3yz) sinnfl; 

r!=0 i'=0 

00 CXD 

"3 {z, 6*) = X! X! (""sni/ cosP^z + sin/3,,z) cosn6'; 

n=0 2^=0 

00 CXD 

M4 {z,e) — (""ill/ sin j3yZ + W4„j^ cos Pvz) cos n^; 



n=0 i/=0 
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oo oo 

U5{z,e) = ^Y^ sin P^z + u-^^^ cos P^z) cos nO (14) 

n=0 u=0 

oo oo 

(r, 0) = 5Z ^-■^ (^) ^) = E (^) 

n=0 n=0 
oo oo 

iy9 (r, 6*) = iy9„ (r) cos n^; (r, 9) — ^I^n (r) sin n9 

n=0 n=0 
oo oo 

Wen (r, 0) = ^ VF-e,„ (r) sin $ (r, 0) = ^ $„ (r) sin n0 

oo oo 

(r, z, 0) = ^ P.^ (r, z) cos 716*; P° (r, z, 61) = ^ P° (r, z) cos n9 

n=Q n=0 

where 

; ' Po' 

In relations H14() "+" is the symmetric and "-" is the non-symmetric components 
of general problem solution. 

Let us exemplify the solution of wave equation ^ which satisfies the Zom- 
merfield's condition in the form of Helmholtz-Huygens integral (6J 

P-(..,0)^£(P-Ug-^UG).<,. (15) 

where dao is element of dissipating surface area ao', 

G — (47ri?*)~^ exp {—iujR*) is the fundamental solution of Helmholtz equa- 
tion; R*^ = r'^ + Tq — 2rro cos (0 — Go) + (z — zq) is the distance between 
point of observation (r, z, O) and point (ro,Zo,6o) on the surface (Tq; is 
the derivative with respect to outside normal of surface (Tq in point (r, z, Q) ; 
(ro, Zq, 0o) is the value of dissipated pressure on surface (Tq- 
Let us arrange the defined representation of function G in cylindrical coor- 
dinates with divided coordinate z 



1 °° r°° p-^\^-^o\ 
G = — V e„, cos [m (0 - 9^)] / A J„, (Aro) J„ (Ar) dA (16) 



where 88= i '^-tj — !p 't' , Jm () is Bessel's function of the first kind; 
1^ zyt^ — A"^, it a; > A 

eo = l;£m = 2;m > 1. 

Let us assume on the surface of closed cylindrical cover the distributed pres- 
sure P*^ (ro, zo, 6o) and its formal derivatives are presented in the following 
form: 
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a) by Fourier's series on the side surface of cylinder 

oo oo 

(zo, do)^Y.Yl ' cos/3,zo + fn^ sm/3,zo) cosnO (17) 

Q OO OO 

—P'^izoJo) = (/n^"" cos /3^,zo + sin /?^.zo) cos n6' 

b) by Fourier-Bessel's series on plates ends zq — ±1 

oo oo 

P"" {zo, 6'o) = X! X! dlfjn hjnro) COS n^o 
n=0 i/=0 

C30 oo 

— (zo, ^o) = E E (7,n'^o) COS nOo (18) 

n=0 !y=0 

where jjn is the zeroes of Bessel's function derivative ( J'j (^jn) = 0) 

Let us write the distributed pressure P'^ as the sum of symmetrical and 
antisymmetric components relative to z = 

P" {r, z, 9) = P^+ (r, z, 0) + P"- (r, z, 6) 

Substituting the expansions l(T7|) . ((TH|) into l(T^ and integrating using 
we get the following expressions for the distributed pressure to region 

a) symmetrical component 

{r, ^) = ^ E f^^ r i^i'^'i^'^ [cos/3.. - (-1)"^ ch (xz) e-'] dA- 
-tE/-^ ^2 [cos/?.. - (-l)-^ Ch {>cz) e-^] dA- (19) 

-T E (7.») ^^4^^^^e-'./. (, >.z ,) dX. 

b) antisymmetric component 

Jo ^{^ +K) 
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__ \ ^ f*e- / 



AJ„ (A) J„ (Ar) 



OC /. 

oo „ 
^ ,._n "'O 



^LA^^l^^^e-'ch{>cz)dX± (20) 

Tin 



e sh (i >r2; i) dA. 



where Eq = 1; = 2; n > 1. 

For the region (r > 0, |z| > /, < 9 < 2tt) expression for (r, z) and P^^ (r, z) 
will look 



- EZo fnZ r ("1)'' (>.0 e-|.|rfA+ 



+ Er=0 fn^Jn (7.n) T ^^^^^^^ (-0 C^^^^dX- 



[ ^EZo9:rJnil^n)jr%^^ch{>cl)e-\^\dX J 



(21) 



+ ET=0 fn.Jn (7,n) ^ ^^^'^ ^^^^ e~-\MX- 

(22) 

+ when z > I 
- when z < —I. 

Functions P^^ {^tz) and P^^ (ro,zo) and introduced in series p7l) - (|18l) un- 
known coefficients fn^,g'^ and functions P^^(ro,2o) are tied by relations 

laiHiini 



1 



fnt ^ ^ j (r, z) \r=i COS {I3yz) dz; eq ^ ^;e^ ^ > 1; 



P,j (r, z) \r=i sin(/3^z)dz; 



(23) 



9n, = ^"J / K [r, z) \z=i Jn hjnr) rdr] e„j = 



27. 



(7|„-"2)-/^(7jn)' 
Substituting the expression p9|l . (|20|) into (|23|l we get algebraic equations 



u=0 



*1±\ 

3=0 



0; 
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where 



^Nvn 



j=0 



j; (A) j„ (A) 







£^/(5^^ (w^ - Pi) 



dX; 



dX- 



G:^•;=--/r^(7,n) 



j: 



^ AV; (A) J„ (A) (-1)^ (1 - e-^-') 



dA 



-dX: 



AV;(A)J„(A) (-irx(l-e-2-0 



<2+ 



Jo 



ijn 
oo \3 Ti2 



aV;^(a)(i 



-2>fi\ 



^^dA-^- 



■2xa 



X^J^ (A) (1 + e- 

°° A^ j;^ (A) (-1)'^ (1 - e-^"') 
(72„-A2) 



dA; 



dA; 



*2+ , ^- AV;(A)J„(A)(-l)-(l-e-^-') _^^ 

- Jn {inn) I — / „ ■„^ -dX\ 



^nvti — Jn iliJ,n) 

Jo 



fL.. = 



Jn (A) Jn (A) 



2 .2^ (-ir^/3./3.AMl-e-2-0 



71*1- 



-I 







(-l)'^+^/3.(l-e-^-') 



dA; 



^nj/n ~ Jn iljn) 



^njn ~ 



oo \2 



PQO ^2 
(7jn) / — 

Jo 



A^j;(A)J„(A) (-If /3^(l-e-^-') 
(7^„-A2)(x2+/J2) 

AV; (A)J„(A) (-If /3^(l-e-^-^) 



rfA; 



dX; 
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AV;(A) (-ir/3.(l 



e 



AV;(A)J„(A) (-ir/3.(l-e-2-') 



oc \ 2 







^2- r . w . W AV;j^ (A) (1 + e-^-') 



/o (7|„-A2)(7^„-A2) 

A3j„2(A)(l-e-^-') 
/o >^(7|.-A2)(72„-A2)' 

is the Kronecker's symbol. 

Let us consider a question of derivatives representation from series H14(l also 
by Fourier and Fourier-Bessel series, as during the solution of equations ©-JHll 
we need to differentiate it. 

Let us assume that even function / (z) and odd function ip (z) and its deriva- 
tives could be expanded into Fourier series in : 

oo oo 

(z) = J2al'' COS P^z;^^'"^ (z) ^ J^^f ^in f3^z; 

oo oo 
J(2k+1) ^ ^„2'=+lsin/?,z;<^(2fe-l) (,) =Y^bf-' COS P,z. (25) 
v=Q i/=l 

CoefHcients al^ ,al^'^ ,bl^'^ are defined by formulas [101 El [HI 



a?,^ = ^ /' [z) cos (i^zdz-hf f ^(2/0) {^z)sin(5^zdz 

2 J~i 2 J_i 



^2k+i ^ /(2fc+i) (z)sin/?,zdz;62fe-i = ^ /' (^(^fc-i) (^)cos/3,zdz (26) 

2 J-l 2 J„; 

Integrating H26|l by parts and taking into account the evenness and oddness of 
derivatives from / (z) and (p (z) we find that 



9k 2£jy 
af ^ — 

2 



2fc+l _ o 2k.i2k _ _n 1.2k-l. 

During the solution of equations we present the function P° (r, z) as the 



series 



E [^„°.+ (^) cos P^z + P„V (^) sin P.z] = P° (r, z) (28) 



i/=0 
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The second derivative from P° (r, z) taking into account 12 5|) . H27|l we write 



as 



|^P° (r, z) = £ ('|-P°+ (r, z) U., ^£;il_i)! „ j ^^^^^^^ 



E 



dz 

2(-ir/3. 

I 



I 



P„°+ (r,z) -/32p°;(r) sin/3,z 



(29) 



Let the dissipated pressure and its derivatives on the interior side of the 
plates are defined as foUows 

Pn^ {r, l)^Y. 9nfJn ilinv) ; (r, = E glfJn hjuV) (30) 

3=0 j=0 

Taking into account p8|l - H3U|) wave equation (jSJ for the symmetrical and anti- 
symmetric components obtains the following form 



92 


1 d 


/ 2 




r dr 


I r2 



(tj"'') ; (31) 



92 1 d 

Qj.2 J. Qj. 



+ -— - I — ->^2 



J=0 



^ 9nj Jn hjnr) ] 



where 88, 



General solutions of H31|) are 

OO 



(32) 



3=0 



where 



Pnu {r) = fny Jn (>^2r) + ^ C„^^-5°j- J„ (7j„r) 



2£.(-l)- _ ^ 2/3, (-1)^ 

"'^J ; /'o^2 ^,2 ^ ' ^ni^J 



From we obtain that pressure (r) and ^P,"* (r) on the interior 
surface of cylindrical cover [z — 1) are 



^^Pn^i^) = /°J(>^2)J;(>^2). 

OO 

p:tw = /„°^^n(><2)-E^-.5:r^«(7,n) 



(33) 
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Let us exemplify P,jj^ from (|32|l by Fourier-Bessel's series and equate it to 
coefficients (|^ . We get 



oo 



o- _ /■0+ (~1) Jra (^2) Jri (Tjn) ^,0+ nr^ , -, ^+ . |'04^ 

~ 2-^ ^,2 _ ,^2 "ly / , \ ^) ^ni^/ji l^^y 

7,n ^2 



iy=0 'J" ^ v=Q 



00 



*o- _ \ " ^0+ (-1)'^ ^2 (>^2) Jn {ijn) ^n^ „o+ ST^ f ^\^n- /'Qt;\ 
3rii — Z^Jnv 2 _ 2 ^rij ^) ^nvii- \^^) 

Let us represent the falling pressure P" (r, 2;, 0) and ints normal derivatives 
on the surface of a cover and plates a the series analogical to dJ, H18|) . H19|l: 



(r, ^, 6*) = ^ (r, z) cos nO; 

00 00 

(1, z) = ^ cos/3,z + sin/3,z) ; P„"±^ (r, Z) = ^ 5,^:^ J„ (7,„r) ; 

^P„" (r, z) 1,^1= ^ (/,:^ cos + /,:r sin/3,z) ; (36) 

Let us present the flat wave as the expansion |Blll5| 

oc 

P" (r, 2, 0) = Po ^ e„i^" J„ (a;sincp*r) exp (— iwz cos cosnd; 
3=0 

and integrating using formulas that define the correspondent series components 
H36|) we obtain 



ra+ (-l)''2e„i "J„ (wsin(^*)wcos(^* 

— COS^ if* 



„a- 2(-l)''£„i "J„ (wsin(p*)/3^ 

/n,. = M 2 — ) ; 

PJ — LJ^ COS-^ if* 



9n^ 5 2 ■ 2 , ^cos(w/cosv5 ); 

7J„ - sm 

S^J = T 2 ■ 2 , sm {ul COS (/9 ) ; 

7j„ - sm (p* 
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j^*a+ --n ■ * Ti / ■ *N 2C0S((5* (-1)'' W . , , *s 

fnu = ^ sm V3 J„ (w sm <^ ) ^2 _ ^2 ^032 y,* ^""^ ('^^ '^^^ -> ' 

/*°" = £„i""w sin (^* j; (w sin ip*) — J" sin {wl cos (^*) ; 

— cos^ (/?* 

= -i^ cos If* g^;j';gZ'' = iuj cos ip* g^^J . 

To find the solution of inhomogeneous equation we define the auxihary 
function 

F± = AoW^± + Ag±; (37) 

^^2 

Taking into account (|37|l to define (r, 0) from (O we get the following 
differential equation: 

(A2 - aiAo - as) (r, 9) = {PI - g± (r, z) ; (38) 

where 

Q Q 

Pi = "2 + "i;p^;P2 = "2-^; 

oc oc 

5^ (r, z)^^Yl (9nf + 9nt - 9nf) Jn {Ijut) COS nO; 

n=0 j=0 

The solution of equation H38() we write as follows: 



i^^ (r,6l) = ^F^ {r)cosne 

ri=0 

CXD 

(r) = Atjn {lir) + A± J„ (72r-) + ^ A3„jff.„, J„ (7jnr) (39) 



/ai± Va?+4a2 , -/3i*7j„ - /32* 

71,2 — y 7i ;^3ni - ^ 2 ; 



^iri'^2ji unknown coefficients. 

Using H39() and H37|l to the normal shifts of the plates we obtain 



7l 72 j^Q tjn 



The solution of equations ©-(jSJ) will be rewritten as follows: 

*n (0 = ^fnJn {hr) ; (41) 

(r) = ^tnJn {kir) (42) 
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^„ (r) = *f„ J„ (fc2r) (43) 

Substituting expansions p4|l into equation taking into account H25|) , H27|l 
we get the algebraic equations relatively to components of shift vector compo- 
nents: 

atunju = bf„ + Sisg^^; af. = a^-; i, j = 1, 2, 3, 4, 5. (44) 



where 



''12 



a± - (1 + a') (-n2 - a^Pl + ^l) - a,(3l - ^; 
a" (n2 - aiPl + xl) -—■af^ = {l + a^) (l + xl) (±n) ; 



= Pu02 {+n) ; afr-, = 0; 0^2 = ( -n^ - (TiPI - ^1 + ^ ) - ^1 ; 



^23 = [(1 + ^1 + a^] (i'^) ; «2i = 0; ajg = a^cr2/3^ (±n) ; 



= - (1 + a ) (Tin - + — ; 045 = a ai - + ^ j ; 



^tn ~ ^^ll'^Ui + ^12'''*2ti + '^M^in! ^2ri = ^21 ""in + ^22^2^1 + ^25'"5n; 
^3n ~ ^33"^^ + ^34"4n + ^SS'^Sri + 9Xw^ ^tn = ^M'^tn + C'45'U|^"„; 
^s'n = ^54 ""In + C'sS^Sn! ^In ^ ^ll^ln + ^12'^2n + C'i4U4„; 
^2ri — ^21^1n + ^22'^2n + C'25M5,i; ^Sn = CsS^Sri + ffrni/; 
~ ^il'^ln + ^43^311 + ^AA^in ^45 ^5n' 5rr!i/ ~ /rw + /ni^ ^ fnt'-> 

+ _ _ 2e.(-l)'' 



^5n — ^b2^2n + ^53'"3ji + C'54U4„ + Qs^Sn; C*!! — 

ri+ _ 2 ( — 1) _|_ ^ 2£y(— 1) ^ _ ri+ 
^12 — '^1 J : ^14 — '^^n J , C-21 — Oi2' 

2£.(-l)'',^+ _ 2, ^2£.(-l)^^+ _ 2 2e.(-l)^ 



6^2 = ~a CT2 ^ ;C^5 = a cr2n ^ ; C J, ^1:4- 

(-l)"^ . ^+ _ ^+ . ^+ _ « 2£.(-l)'^^ _ „2. 2e. (-1)^ 



C34 — -^'01 ; ; C'ss — (^33; (^44 — /^i- ; ^ Q5 — a (3v- 



r+ — r+ ^ r+ ■ r- — ^ 2(-l) _ „2 2(~1) „- _ „ 

'-'54 — "-^55 + ^45! '-'ll — '^'^ ^ ' ^12 — ^ ' ^21 " ^12; 



-f- _ ,,2 2( 1) . _ ^ ^ 2( 1) _ — . 

-■33 — ^ 1 '-'33 ^ ~cr2" ^ ,'-'53 — ^33,054 — O55 — O45, 
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C52 — «^C'4i; Cf^ — ±C^; 



From 1031) we define f/j^j,, as: 



(45) 

where 

-^1^1/ ^ ('^n^lj' + ^^21 ^2;-) ; -^2niy = '^T ('^12^1j- + ^22'^2j) ' 



1 A + 

]+ _ _!_ /'/^+ 4+ ^ . r J+ _ _Jl. 



^2nv — 12"^li + ^22^2j + 52^5jO 

5 

^3jii^ = (C 33 ^3j + ^43^4] + ^ 53^5j] 

L\nv = "T^ (C 44^4j + ^54^5j) ! -^W ~ (^^ 45^4j + C'sS^Sj) 

^5 ^5 

Lilu = = det 14 1 ; J = (1, 2, 3, 4, 5) ; 



Af- are algebraic adjuncts of elements a^. 

To improve the convergence of the obtained solution H45I) we use the asymp- 
totic properties of expansion coefficients of Fourier series of the desired functions. 
Basing on (|25|l . (|27|l we could define the following relations to the even functions 
/ (z) and odd functions (z) and its derivatives represented as Fourier series 
on [—1, 1] , 

a) for even function / (z) 



f {z)=^f u COS (3yZ+ ^ 



2£. (-!)'• 



v=N+l 



I 



f'ii) /"'(0_^r(o 



f^L^ 1^1' f^V 



COS f3iyz; 
(46) 
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N 

E 



cos I3,jz+ ^ 



2e.(-ir 



u=N+l 



I 



cos /3^z 



TV 

E 

j/=0 



(/3^/' (0-/"'(0) 



COS/3iy2: + 



+ E 



2£. (-1)^ 



v=N+l 



I 



PI Pt 



etc. 



b) for odd function 



N 



(z) = ^(Pi, sin /3^z; 



iy=0 



jz) 
dz 



N r 



E 



2(-ir 



^(0 



cos p^z+ J 

i^=N+l 



Pt 



COS PuZ 



(47) 



E 



d'^ip [z) 
dz^ 

oo 

- E 

i/=JV+l 

(z) 

2 (-If 
I 



N 

E 

i/=i 

2(-ir 



(/) (/) 



PI Pt 



cos P^z : 



d^(p{z) 



N 

E 



(,3»^(0-/W (0) 



C0S/3,yZ + 



cosPi,z+ 



V" (0 (0 



u=N+l 

Let us introduce the symbolism 

Xi — u'lj^] X2 — u'^i, ~ Ug^; X^ = u^^; X^ — u^^; 



dz^ 



^1 ="l«;^2 ="2n;^3 =W3n;^4 = ; ^5 ="5n; (48) 

Substituting obtained solutions of differential equations ©-I© in^to the bound- 
ary conditions we get the system of equations relative to unknown con- 



stants (|48|l with asymptotic H46|l . (|47|l for Uj^^ taken into account: 



(49) 
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where 



N 



N 



N 



N 



< = E (-1)^ ; "^2 = E (-1)'^ ; ^^3 = E ^'^n. i-ir+Bi;at, = E LlL (-1) 

i/=0 i/=0 1^=0 1^=0 

^ 21 °° 1 

(^i5 = Y,^5Li-'^T''at5 = l-,afg = ~kiJ!^{ki)-,af^Q = -nJn{k2)-,B* = - ^ 



i/=0 

N 



N 



N 



!y=JV+l 



4i = E (-1)^+^1*; 42 = E (-1)^ ; 43 = E (-i)'' ; 44 = E ^4+. (-1)' 



i/=0 



i/=0 



^^26 ~ 



N 



Jn (71). 
71 



"^^ ; 4io = ^2 (fe) ; afg = -nJ„ (/ci) ; 



72 



TV 



N 



N 



45 = E (-1)"^ ; 4i = E (-1)'' ; 42 = E (-i)''+sr; 43 = E ^3+. (-1)' 



v=0 

N 



v=0 

N 



i/=0 



v=0 



L/=0 



i/=0 



44 = E (-1)'' ; 45 = E ^3^>/ (-1)'' ; 44 = i; 49 = (^i-^n (''^i) - Jn {kl)) ; 



nJn {k'2) k2Jn (fc2) 



h 571 j; (71) ; a 

71 



— + ^72 J;(72); 
72 



Ro 



N 



E4t(-ir+i?i* 



.!^=0 



i?0 



AT 



E^?-(-ir 



1^=0 



TV 



"58 



N 



X+3 = E (-1)"^ + E ^3t (-1)^ 



TV 



Ro 



u=0 



i/=0 



Ro 



N 



u=0 



N 



.v=0 



N 



44 = ^1 E4t(-ir/3. 



N r 



45 = ^1 E 



2N 



N ^1 

t^o ^0 ;^ J 



i5t(-ir/3.+^ 



+ ^-E^5t(-ir+^E^it(-ir 



TV 



i/=0 



-Ro 



i/=0 



= -A'57i j; (71) ; = -A'572 j; (72) ; fles = -'^■^n (^3) ; 44 = 45 = ^'i 



^710 



{i-iy^2)R<y 

-£/2/l2 



Vo2f^ Jn (^1) I i 

[nfc2 (^2) + i^02nJ„ {k2) - i^02nJn (fc2)] ; 



«81 - -°1 2 ' 85 - 2 R^' 
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E2h2 



(1 - i/q's) ^0 

-E2/12 r .2 



[-nkiJ'n (fci) - nkiJ'n (fci) + nJ„ (fci)] 



(1-^^02)^0 

' N 



[-n^Jn (fc2) - fci-^n (fc2) + ^2 j; (fcl)] ; 



i?0 



.i/=0 



JV 



a+3 = ^1 ^ E ^st (-1)^ /3. + ^ E ^3+. (-l)"^ ^ ; 



N 



i/=0 



AT r 



44 = £'l E 



2Ar 



JV 



i/=0 



JV 



JV 



i/=0 



a^6 = -Di 



0^7 = Di 



45 = ^i<| E (- + ^ E (-i)"^ I 

(7?5 - 1) j;; (71) + 1^02 ( 571 - ^^'j j; (71) + ^^02n2 f^-ii 



{^IS - 1) (72) + Z^02 ( ^72 - -\ J'n (72) + (^-^ 

72/ V 72 



48 = --^2 j ^3 Jn (A:3) ; afoi = ^2 ^ 2R0 J ' ""'^ 

4o4 = --Di 4o6 = -'D2n (1 - 1^02) (^Sji - j; (71) ; 

4o7 = -D^n (1 - ^'02) (^^72 - j; (72) ; a^3 = 1; 



1 - 1^02 



JV 



JV 



JV 



JV 



«21 = E-^l»>'(~^)''5"22 = E^2n>/(-l)'';a23 = E^3ni.(-l)'';a24 = E ^4n^ ("l)''; 

iy=l l/=l l/=l l/=l 

JV JV JV 

= E ^5n>. %! = E ^l"-^ ; «52 = E ^^n:. ' '^Sl = «l2=l; 

f=l f=l I/ = l 

TV TV W 

«51 = E ^1- (-1)' ; «52 = E ^tn. i-^r ; «53 = E ^3^. ("l)'' I «54 = E 



L/=l 



L/=l 



1^=1 
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N 



^5n^ (-1) S^ei = -^i^R^'^ea = ^i;r^'"65 = 



^L?„,(-ir/3. + ^Lt,(-l)'^ ; 

i/=l i/=0 J 

{AT N ^ 

^ L^-, (-ir /3. + ^ ^ L^-, (-If I ; 



-^73 



= Ai (E ^L". (-1)^ + ^ + E ^L". (-l)-- 1 



TV 



TV 



%4 = Ai E^4;.(-ir/3. 

U=i 

{AT Af 
EiL.(-ir/3.+E^w(-ir ; 

afe = -Ai |e (-1)'' + (-l)T 5™. I 



1 - 1^02 J n 



1/=1 



.i/=0 



2Af 

T 



«82 = Bl 



«83 = -Si 



2 I Ro 



1-1^02 \ n 



f^L^-,(-ir + f^Li-,(-ir/3J; 



I/=l 



i/=0 



I -1^02 \ n 



«85 = -^1 " 



1 - 1/02 I n 



"lOl = -D; 



1 - J/02 I 1 



<^in2 - -^1 



2 1 

1 - 2^02 



{ E ^in. i-^y + E ^3;;. (-1)'' /3. 1 ; 

{AT AT 
-^E^4;.(-ir+E^i-(-ir/?4; 
1^=1 1^=1 J 

■ AT 1 ^ 1 

.i/=l J v=0 J 

Af AT 



Ro 

N 



"92 = ^1-5-^' "94 = -'^i; 

-Ko 
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1 - 1^02 

2 

1 - 1^02 



{N 1 ^ ^1 



^104 - ^1 — 2 



n 
Ro 



N 



N N >, 

^ L^;, (-ir - — ^ Li;, (-If 



'-105 



1 - 1^02 



{N 1 ^ ^ ^ 



N 



N 



A 



3nj 



— E/ ^'inufrrii' + E] ^inudtnj'^ * — 1: 2, 3 10 

4 SL 



•/n iljn) ; Cgn, — 



^5+ , Vain ^1+ vm 3+ 
/to /to 



J-4+ , J- 2+ 



Ro 



i:'2 



t^n {ijn) ] 



A _ .,2 Q A s;_ 

A3n, 7,-„^A3„,- ^^^^^ 
Tin 



1 - f02 



Solution (gni will be 



AT 



— E ■^nvjf^nu + E ^nu]9fnp (* ^ 1' 2, 3. ...10) 



(50) 



where 



i/=0 



7± = 



J=0 



E 



,± 



k=\ ^10 



For fJsniy and U zny. we obtain 
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To improve the convergence of the obtained solution we will single-out the 
characteristic features in unknown functions. It is known that ^J^] approach- 
ing to angle circles on surface r = 1 and z = I expression for the pressure P^^ 
and its derivatives must have the following singularities 

a) for the pressure 



(zo) = \]ip~zlf- Pr (zo) = Bn^o ^J{l^-zl): (^o) = Dt {I - r^f : 

(52) 

b) for the derivatives from pressure 

9P„^+ (zor) _ Al OPr (zpr) P^o 



=1 



Expanding H52() . (|53|l into Fourier's, Fourier-Bessel's series corresponding to 
(O we obtain 12 , 15 



An Y - zlf = ^ /,"+ cos /3^Zo; Sn2:0 \/ ('^ - ^o) = ■^"'^ '^^'^ ^''^O' 

i/=0 i>=0 



I °" /I* °° 

u=0 ^0 ^=0 



where 



.+ ^ AjV6r(5/3)r(i/2)r,/6 {(3j)2'/\ _ i?;±£„,r(2/3)r„+7/6 (7,„). 



.± ^ At£n,2^/'r(5/3)r„+5/3 (7,-„) ^ A„r(2/3)r(l/2)Ji/6 {pj)2^i^ 

^ (5/3) r (1/2) Ji3/6 2^/6^/6 ^ A„r(2/3)r(l/2)Ji/6 {PJ)2^'^I3J 

^5/6^7/6 ^7/6 

Unknown coefficients /^^, c?;;^, at large v,j{v,j > N) behave 

itself like coefficients of known expansions (|54(l . i.e. we could write 

j'eit p+e± n^"*" n**^"'" 
rei J riN ts±. r^ezL nN j-^si, ^e± JriN_ si , /^^^^i ynN_ ^s:L (^^\ 

/nAf JnN 9nN 9nN 
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Thus, if in the unhmited systems of algebraic Hnear equations H24|) unknown 
coefficients are replaced with formulas H55|) than elements of iV-th columns in 
matrixes will be: 

_(-!)" /•/ aV;(a)j„(a) . ^ 




j=N+l 



" <r i 1 ^ ^ ^^^^ ^ 

p,2+ Jn(7,») /-J AV;(A)J„(A) , (-ir/*g+ l 



^„^7V - p+ / \ 9 _X2 ll ^ e ) ^2, «2 f ^'A, 



(56) 



^.2+ _ Jn (7,n) /• J AV;^ (A) E^TM^UX 



^2+ _ Jn {l,n) I I A3 (A) . n^,. f^. Q^pn {ijn) I 




A-^K (A) / _ -2>,l\ Sr^ j-^y Pvfni 



■d\: 



AJ^ (A) Jn (A) , ^ -2>cl\ (-1)'" Pvfnl 1 

AV„(A) J„(A) _2^;. ^^TjMTjrOl 

i_ _ (-1)^/3^ /■ J AV;(A)J„(A) ,^ 



24 



^.2- Jn iljn) f \ (A) J„ (A) ^ {-^rMZ- \,.. 

/:^- yi-(7^-A^) ..^+/?? /'^' 

^.2- _ (7.») /• / (A) ^TMZ^Lv 

- / \ —7-2 ToT 11 - e ) 2^ 2 _x2 f o'A, 

^*2- _ J» (7j») /• J AV;^ (A) , ^ gn7-/«(7jn) [ 

<^n,.iV - -i— / \ ( 2 _^^2) + ^ ) ^2 _;^2 f "'A. 

5„Ar Ujn j 3=N+1 ^J" J 

where 

2i/6^i/6r (2/3) r (1/2) j,/6 (^0 (-1)'^ ^ f> fnt-'i-ir 

2w/6r(5/3)r(i/2)/7/6H) ^ v-+(-ir _ V 

i<2 I /92 



2^/3AV3r (2/3) J„+5/3 (A) g^lf J„ (7,-„) _ J»(7,») 

(A2-„2)J2(A)J„(A) ^4^^ 7|„-A2 7|„-A2 ' 

22/3^-2/3r(5/3)r(l/2)J„+2/3(A) ^%;;±J„(7,„) _ ^ J„ (7, 



(A2-n2)j2(A)j„(A) A. ^2^„A2 -^A.^^ ^2^_A2 ' 

2 (2/3) r (1/2) [>.^/;^, (>./) - i/e/ye {^i)] f*s- (-i)- _ - /--(^if/?. 

27/6^i3/6r (5/3) r (1/2) [,,^4^^ (,,/) - 7/6/,/a /^+ („i)- _ - /^t+l-lf/?. 

Substituting the obtained values into condition of hydro-resihent contact (0 
and taking into account (|24|l , 134(1 , ((35|l , 1)56(1 we obtain such unhmited systems 
written in matrix form to evaluate the unknown coefficients 

re± r*e± e± *e± ro± r*o± o± „*o±.l 0^012 N)- 



NxN IxAf NxN IxW iVxJV IxAf iVxJV IxJV 

(57) 
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NxN IxJV JVxAf IxW NxN IxN NxN IxN 

(58) 

ATxiV IxJV NxN IxN IxW IxTV 



ATxJV Ix^^ WxAf IxN IxN IxN 



{Rt)Jftn}, + {RL),,{9tn}, = «),,{/r^}.; (61) 

NxN IxAf NxN IxN NxN IxiV 



NxN IxAf ATxJV IxAf NxN IxW 



where 



NxN IxN 



TVxAT IxTV TVxAf IxJV IxW 



jVxAT IxAT 



NxN IxN 



IxN 



R 



1,2,3,4. 



ca:L 

f 711/ 



9l1 



E, 



4:71^11/ 



J nu ' J 71 
■ TV 



9lt 



2N 



N 



X2Jn (^2) . 



E 



4:71 flU 



M. 



^ni^U^fJ-n (-1)*" - 

(-1)" >tf2>/„ Jn {ijn) 



(-1)" K2Jn {^2) Jn {ijn) SfmPu 



lln - A 



(63) 
(64) 



Thus we obtain the closed finite system of 8 (2iV + 1) hnear algebraic equa- 
tions for the symmetrical and anti-symmetrical components of the solution of 
given problem. From the system H57|l - H64|) we define the unknown coefficients 

ve± r*e± e± *e± ro± r*o± o± *o± j = I 2 N 
J nu 7 J nv 7 ynj ^^nj ^ J nu ^ J nu Tifnjiifnj i ^ ^ J 
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Let us consider some particular cases, which are directly obtained from (jSTJ: 

1) If liquid filler is absent inside the cylindrical cover than we obtain the 
soluted system of first four matrix equations (|S7|) with the following mod- 
ifications 

2) For absolute solid (soft) cylinder the soluted system will be formed of 
matrix equations l(57|) . ifsSI) and such two equations 

{/n }i/~{/n {_9n }j~{5n },y~{/n }yi {_9n }j ~ {.9n }j i 

Nxl Ny.1 ATxl Afxl Nxl Nxl atxI Wxl 

(65) 

3) If the cylindrical cover is replaced with liquid cylinder we obtain the so- 
luted system containing H57|) . H58|l . H()3|l . H64|l matrix equations and such matrix 
equations: 



Nxl Nxl NxN uu Nxl Nxl NxN Nxl 



(66) 



Afxl Wxl JVxl Afxl ATxl TVxl 

4) If point pressure source 

is placed inside the cylindrical cover than the system of matrix equations H57|l 
- (|64|l must be changed as follows 

/r+ = /:r = g^J = gT = «; (67) 

a+ _ J^O£t/£nj (-1)" J71 hjn) COS P^Zp _ PpP^enj j-l)" Jn jljn) siu/j^Zp 



PgSnjJl ( 7jn) sin/3t,Zo a+ _ fo£iy£n3 Jn (7j«) COS/g^Zp 



j=0 v=a 



Thus, in this work a new problem solution methodology of contact interac- 
tion of acoustic medium with resilient finite bodies of cylindrical form, based on 
application of boundary integral equations method in conjunction with series 
method with later use of series convergence improvement taking into account 
the particularities of the determined functions. 
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